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ABSTRACT 

We present an X-ray study of the low-luminosity active galactic nucleus (AGN) in NGC 4258 
using data from Suzaku, XMM-Newton, and the Swift/HAT survey. We find that signatures of X- 
ray reprocessing by cold gas are very weak in the spectrum of this Seyfert-2 galaxy; a weak, narrow 
fluorescent-Ka emission line of cold iron is robustly detected in both the Suzaku and XMM-Newton 
spectra but at a level much below that of most other Seyfert-2 galaxies. We conclude that the 
circumnuclear environment of this AGN is very "clean" and lacks the Compton-thick obscuring torus 
of unified Seyfert schemes. From the narrowness of the iron line, together with evidence for line fiux 
variability between the Suzaku and XMM-Newton observations, we constrain the line emitting region 
to be between 3 x lO^r^ and 4 x lO^r^ from the black hole. We show that the observed properties 
of the iron line can be explained if the line originates from the surface layers of a warped accretion 
disk. In particular, we present explicit calculations of the expected iron line from a disk warped by 
Lens-Thirring precession from a misaligned central black hole. Finally, the Suzaku data reveal clear 
evidence for large amplitude 2-10 keV variability on timescales of 50ksec as well as smaller amplitude 
fiares on timescales as short as 5-10 ksec. If associated with accretion disk processes, such rapid 
variability requires an origin in the innermost regions of the disk (r w lO^g or less). Analysis of the 
difference spectrum between a high- and low-fiux state suggests that the variable component of the 
X-ray emission is steeper and more absorbed than the average AGN emission, suggesting that the 
primary X-ray source and absorbing screen have spatial structure on comparable scales. We note the 
remarkable similarity between the circumnuclear environment of NGC 4258 and another well-studied 
low-luminosity AGN, M81*. 

Subject headings: black hole physics, galaxies: Seyfert, galaxies: individual (NGC 4258), X-ray: galax- 
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1. INTRODUCTION 

The low-luminosity active galactic nucleus (LLAGN) 
NGC 4258 is an essential object in our quest to un- 
derstand the astrophysics of extragalactic supermas- 
sive black holes. Very Large Array (VLA) and Very 
Long Baseline Interferometry (VLBI) observations have 
found a set of water masers that trace a nearly edge- 
on, geometrically-thin gas disk ^ 0.2 — 0.3 pc from the 
central black hole (Miyoshi et al. 1995; Herrnstein et 
al. 1999). The near-perfect Keplerian velocity curve 
of these water masers provides one of the strongest and 
most robust pieces of evidence for the existence of extra- 
galactic supermassive black holes (Miyoshi et al. 1995). 
Furthermore, the dynamics of this disk allows precise 
measurements of the central black hole mass, (outer) ac- 
cretion disk inclination and warping, and its distance 
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{M = 3.9 ± 0.3 X 10^ Mq, Z) = 7.2 ± 0.5 Mpc, yield- 
ing rg EE GM/c'^ = 5.8 x 10^'^ cm, Ipc = 5.4 x 10^ rg, 
and 1" — 35 pc; Miyoshi et al. 1995; Herrnstein et al. 
1999). This is the only AGN for which the black hole 
mass, distance, and (outer) accretion disk geometry are 
so accurately known. 

Due to the precision with which the black hole mass 
and distance are known, NGC 4258 can be used as a ba- 
sic test-bed for our models of black hole accretion. The 
overall luminosity of the AGN is small compared with the 
Eddington luminosity of the black hole, L ^ lO^^'LEdd, 
the cause of which has remained controversial. Is the 
small luminosity simply due to a very small mass accre- 
tion rate through a radiatively-efficient disk, as suggested 
by modeling the physics of the maser production Neufeld 
& Maloney (1995)? Or does the disk make a transition to 
a radiatively-inefficient accretion flow (RIAF) at some ra- 
dius as suggested by the modeling of Lasota et al. (1996) 
and Gammie, Narayan & Blandford (1999)? What is the 
role of the jet, and how much of the radiative luminosity 
of the AGN is actually due to the jet (Yuan et al. 2002). 
In the bigger picture, what is the fundamental differ- 
ence among the accretion flows in the most undcrlumi- 
nous galactic nuclei (e.g., Sgr A*), LLAGN and powerful 
AGN? 

Sensitive X-ray observations provide a powerful means 
of probing both large scale and small scale structures 
within NGC 4258 and hence addressing these questions. 
Soft X-ray thermal emission from hot gas associated 
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with the weU known hehcaUy twisted jets (the anomalous 
arms) has been known since the Einstein days (Fabbiano 
et al. 1992; also see ROS AT work of Pietsch et al. 1994; 
Cecil, Wilson & De Free 1995; Vogler & Fietsch 1999). 
However, power-law X-ray emission from the (absorbed) 
central engine of the AGN itself was not seen until the ad- 
vent of ASCA (Makishima et al. 1994; Reynolds, Nowak 
k Maloney 2000). ASCA clearly revealed variability of 
both the absorbing column and hard X-ray flux on the 
timescale of years (Reynolds et al. 2000; hereafter ROD), 
a result supported by short observations of NGC 4258 
by XMM-Newton and Chandra (Fietsch & Read 2002; 
Young et al. 2004; Fruscione et al. 2005). The most sen- 
sitive hard X-ray (> 10 keV) study of NGC 4258 was con- 
ducted by BeppoSAX (200ksec; Fiore et al. 2001) which 
detected the AGN emission out to beyond 50 keV (Fiore 
et al. 2001). BeppoSAX also revealed day-timescale 
variability of the hard X-ray power-law, setting a firm 
upper limit of 250 (5 x lO^'^pc) to the X-ray emission 
region size. 

In this Paper, we present results from new Suzaku and 
XMM-Newton observations of NGC 4258 which, when 
supplemented with survey data from the Swift Burst 
Alert Telescope (BAT), gives us an unprecedented view 
of this AGN from 0.3 keV upto 140 keV. These data sug- 
gest a circumnuclear environment that is remarkably 
"clean" compared with other Seyfert 2 nuclei. Suzaku 
also reveals rapid variability of the AGN emission, allow- 
ing us to set new constraints on the size/compactness of 
the X-ray source. The plan of this paper is as follows. 
Section 2 briefly discusses the data that we utilize and 
the basic reduction steps. We present our analysis of 
the spectrum, as well as spectral variability, in Section 3. 
Section 4 discusses the implications of these results for 
our understanding of the structure of this AGN and the 
origin of the X-ray emission. In particular, we argue that 
the iron line has all of the properties expected if it were 
to originate from the surface layers of the (warped) ac- 
cretion disk. Throughout this paper we quote error bars 
at the 90% confidence level for one interesting parameter. 
All error bars on figures are displayed at the la level. 

2. OBSERVATIONS AND DATA REDUCTION 

In this section, we discuss the new observations pre- 
sented in this paper and the subsequent data reduction. 

Suzaku observed NGC 4258 for a total of 186 ksec start- 
ing lO-Jun-2006 as part of the Cycle-1 Guest Observer 
Frogram (FI: Itoh, US Co-FI: Reynolds). AU four X-ray 
Imaging Spectrometers (XIS 0-3) as well as the Hard 
X-ray Detector (HXD) were operational and collecting 
data, and NGC 4258 was placed at the "nominal HXD" 
aimpoint. Reduction started from the cleaned Version- 
2 data products, and data were further reduced using 
FTOOLS version 6.4 according to the standard proce- 
dure outlined in the "Suzaku Data Reduction (ABC) 
Guide". The standard filtering resulted in 88.5 ksec of 
"good" XIS data. Spectra and lightcurves were extracted 
from all XISs using a circular region of radius 3.25 arcmin 
centered on NGC 4258. Background spectra were ob- 
tained from rectangular source free regions (avoiding the 
calibration sources) around the chip edges. Response 
matrices and effective area curves were generated using 
the xisrmf gen and xissimarf gen tools, respectively, us- 
ing the recommended 400,000 photons per energy bin 



during the construction of the effective area files. We 
also utilize HXD data in this paper. Standard filtering 
resulted in 97.5 ksec of "good" Suzaku-HXD/PIN data 
from which a spectrum was constructed. A FIN back- 
ground spectrum was produced that included the Cosmic 
X-ray Background (CXB) plus the latest model of the de- 
tector background. We do not consider HXD/GSO data 
in this paper due to the fact that the high background of 
this detector makes a detection of NGC 4258 impossible. 

The XMM-Newton data presented in this paper re- 
sult from a continuous 65 ksec exposure started on 17- 
Nov-2006, 160 days after the Suzaku observation. All 
of the European Fhoton Imaging Cameras (EFIC) were 
operated in PrimeFullWindow mode, and the data were 
cleaned using the Science Analysis System (SAS) ver- 
sion 7.1 following the standard procedure outlined in 
the "User's Guide to XMM-Newton SAS". We reject 
data during three background flares (FN count rate 
>50cts^^). The final "good" exposure time for the 
EFIC detectors is 55 ksec. EFIC spectra were extracted 
using a circular extraction region of radius 30 arcsec cen- 
tered on the bright nucleus of NGC 4258, and back- 
ground spectra were extracted using a nearby source free 
circular region of radius 1.5 arcmin. 

The Swift/BAT is a wide field (2 steradians) coded 
aperture hard X-ray instrument which, during nor- 
mal operations, surveys 60% of the sky each day at 
<20milliCrab sensitivity. The BAT spectrum used here 
was prepared as part of the 22-month extension to the 
BAT-survey. The BAT survey spectra are derived from 
an independent all sky mosaic map in each of nine energy 
bins averaged over 22 months of data beginning on Dec 
5 2004. The energy bin edges are 14, 20, 24, 35, 50, 75, 
100, 150, 195 keV. The nature of the coded-mask anal- 
ysis naturally results in a background-subtracted source 
spectrum. As discussed in Tueller et al. (2008), fitting of 
the BAT data was performed using a diagonal response 
matrix which correctly accounts for instrumental system- 
atics in sources with spectral indices similar to the Crab 
(photon index F '--^ 2). See Tueller et al. (2008) for more 
details. 

All spectra were binned to a minimum of 20 counts per 
bin to facilitate fitting. All spectral analysis presented 
here is performed using XSFECvll.3.2. 

3. RESULTS 

In this section, we describe the spectral and temporal 
properties of the X-ray emission from NGC 4258. Given 
its superior signal-to-noise and energy band our discus- 
sion focuses on the Suzaku dataset, although we compare 
and contrast with the new XMM-Newton dataset where 
ever appropriate. In order to extend the energy reach of 
our time-averaged spectral study fi )3.ip . we supplement 
the Suzaku data with Swift/HAT data. Finally, in order 
to study long-term variations in the all-important fluo- 
rescent iron line, we include data from our new XMM- 
Newton observation as well as archival XMM-Newton and 
ASCA data. 

3.1. The 0.3-140 keV spectrum 

The combination of the Suzaku/HlSs, the Suzaku/PIN 
and the Swift/BAT allows us to form the spectrum of 
NGC 4258 over the energy range 0.3-140 keV (see Fig.[T] 
which, for clarity, only shows one of the four XISs). This 
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Fig. 1.— The 0.3-140 keV spectrum of NGC 4258. Shown here 
are data from Suzaku/XISl (0.3-10 keV; black), Suzaku/Pm (14- 
30keV; orange), and the Swift/BAT survey (15-140 keV; gold). 
The spectral model consists of two optically-thin thermal plasma 
component (magenta; unabsorbed apart from the effects of the cold 
Galactic absorption column of A^cal = 1-45 X 10'^'^ cm~'^), an ab- 
sorbed power-law (with intrinsic absorption Nh = 9.2 X lO^'^ cm~^ 
and photon index F = 1.77), and a narrow 6.4 keV iron fluorescence 
line. 

spectrum extends to significantly higher energy than any 
previous study of NGC 4258. It must be cautioned, how- 
ever, that the Swift/BAT data are cohected over a 22- 
month period compared with the Suzaku "snapshot" . In 
this study, we make the assumption that the form/shape 
of the high-energy (> lOkeV) spectrum does not change 
with time even if its normahzation does change. This 
allows us to perform joint spectral fits of the Suzaku and 
Swift/BAT data in order to study the nature of the X-ray 
source and the circumnuclear environment. 

Guided by previous studies of this source, we model 
this spectrum as the superposition of optically-thin ther- 
mal plasma emission with temperature T (described by 
the XSPEC model vmekal; Mewe et al. 1985; Kaas- 
tra 1992; Liedahl et al. 1995), an absorbed power- 
law (with photon index F absorbed by a cold column 
iVu), and an additional continuum component required 
to correctly describe the inflection point in the spec- 
trum around 2 keV where thermal plasma emission and 
absorbed power-law swap dominance. There are sev- 
eral possible identifications of this additional continuum 
component including (1) AGN emission that has scat- 
tered around the absorbing matter, (2) AGN emission 
that has leaked through a patchy absorber, (3) hard 
X-ray emission associated with X-ray binaries in the 
galaxy, or (4) thermal emission from very hot gas as- 
sociated with star formation or the interaction of the 
AGN jet with the galactic disk. Both the intrinsic ab- 
sorption and the Galactic absorption (from a column 
Nu = 1.45 X 102°cm"2. Dickey & Lockman 1990) were 
described using the phabs model. To obtain a good de- 
scription of the soft X-ray data, we require that the el- 



emental abundances are allowed to vary relative to so- 
lar values (as defined by Anders & Grevesse 1989) in 
two groups. Group A (with abundance Za) consisting 
of {C,N,0,Ne,Na,Mg,Al,Si,S,Ar,Ca}, and group B (with 
abundance Zb) containing {Fe,Ni}. 

In our Canonical Spectral Model, the additional con- 
tinuum component is modeled by a power-law component 
with photon index equal to that of the main (absorbed) 
AGN powerlaw = F) that is unaffected by the intrin- 
sic absorption. This is an appropriate spectral model 
to describe the scattering or leaky absorber scenario. 
This model provides a decent description of the spec- 
trum (x^/dof = 4506/4003) with best-fitting parameters 
F = 1.75l°:°t = {9-2toi) X 1022 cm-2, kT = 0.54 ± 
0.01 keV, Za = O.49t°;J°Z0, Zb = O.27l°:°^Z0. The 
normalization of the additional continuum component 
relative to the main AGN powerlaw is / = 6.0 ± 0.4%; 
this can be interpreted as the scattering or leakage frac- 
tion. If we allow the photon index of the additional con- 
tinuum component to be free, the fit does not improve 
and we find that r2 is very poorly constrained. Using 
this spectral model, we deduce an observed 0.5-10 keV 
flux i^os'Lio = 9.7 X 10~^2 ergcm~2 and observed 

0.5-10 keV luminosity L^^i_io = 6.0 x lO^^ergs"! . Re- 
moving the effects of absorption within the model im- 
plies an intrinsic 0.5-10 keV luminosity of ijfs-io — 
1.4 X 10"'^ ergs^^ . The best- fitting normalization of the 
BAT model is 48% that of the Suzaku model. We inter- 
pret this as true variability, i.e., the Suzaku observation 
caught the source at a time when the high-energy emis- 
sion has twice the normalization of the 22-month average 
spectrum. 

The second continuum component can also be de- 
scribed by an additional thermal plasma component; 
however, while a good fit can be obtained (x2/dof = 
4501/4004), a rather hot {kT > 5.5 keV) and low- 
metallicity {Z < 0.26^©) plasma is required. Thus, since 
it is not clear from where such a plasma component would 
originate, we prefer the power-law description of this ad- 
ditional continuum component. 

Our Canonical Spectral Model fit leaves a line-like 
residue in the 0.5-0.6 keV range that likely signals soft 
X-ray complexity beyond the simple one-temperature 
thermal plasma model. Adding a second plasma com- 
ponent with an identical abundance pattern but differ- 
ent temperature {kT = 0.22 ± 0.02 keV) leads to a sig- 
nificant improvement in the goodness of fit {x^/dof — 
4360/4001). However, a robust exploration of these 
multi-temperature solutions is not possible at XIS- 
resolutions (e.g., the abundances are poorly constrained) 
and we defer further discussion of this to a future pub- 
lications in which high resolution spectral data from the 
XMM-Newton/RGS and a long Chandra/RETG obser- 
vation) will be discussed. 

The combination of the XIS, PIN and BAT data allows 
us to examine the hard X-ray spectrum of this source in 
more detail than previously possible. For the hard-band 
study of this paragraph, the data below 3 keV were for- 
mally excluded from the analysis to prevent their high- 
statistical weight from biasing the high-energy fit. No 
significant deviations from a pure absorbed power-law 
are detected above 3keV. Of course, in order for the 
observed powerlaw (with F < 2) not to possess a di- 
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vergent energy, it must cut-off or rofl over at some high 
energy, ff the power-law has an exponential cutoff, the 
characteristic e-folding energy is constrained by our high- 
energy data to be iSfoid > 124 keV. If we instead assume 
a pure-power law with cold X-ray reflection (modeled us- 
ing the pexrav code; Magdziarz & Zdziarski 1995), the 
constraints on the "reflection fraction" TZ are very depen- 
dent upon the assumed inclination of the reflector. If the 
reflector has a slab geometry with a very high inclination 
(i > 80°), as expected if we identify it with the inner 
accretion disk, then these data provide no meaningful 
constraints on the reflection fraction. For more face-on 
reflection (e.g., the surfaces of discrete cold clouds or the 
inner wall of a cold torus on the far side of the X-ray 
source), we find TZ < 0.43. If we allow both reflection 
and an exponential cutoff simultaneously, the limits on 
the e-folding energy become weaker (i^foid > 67 keV) but 
constraints on reflection are essentially unaffected. The 
implications of this result for the circumnuclear environ- 
ment in this Seyfert nucleus is discussed in §4. 

An identical analysis of the 0.7-10 keV XMM- 
Newton/EFIC (PN and MOS) data (supplemented with 
the Swift/JiAT spectrum) gives a similar picture, al- 
though there are some quantitative differences. Firstly, 
the flux of the thermal plasma emission is lower by a fac- 
tor of two, a simple consequence of the fact that much of 
this emission lies in the anomalous arms and is outside 
of the our EPIC extraction region. The temperature and 
iron abundance of this emission {kT — 0.58 ± 0.01 keV, 
Zb — O.17lg'Q5^0) are very similar to that derived from 
Suzaku, although the light metal abundance is slightly 
lower {Za = 0.20lg j!)gZQ). The most robust differences 
are with the parameters describing the absorbed power- 
law; the power-law is flatter (F — 1.65 ± 0.07) and less 
absorbed [Nh = (7.7 ± 0.5) x 10^^ cm"^] in the new 
XMM-Newton data as compared with the earlier Suzaku 
observation. Accompanying these changes is a decrease 
in intrinsic (unabsorbed) 0.5-10 keV luminosity from 
= 1.4xl04iergs-i to L^g.ig = GxlO^Oergs-i . 
We will see below that these long-term changes are in 
the same sense as short term variability seen within the 
Suzaku observation and may be revealing aspects of the 
spatial structure of the X-ray source and absorber. 

3.2. The Suzaku detection of a weak iron line 

The width, strength, and variability of the 6.4 keV 
K-shell fluorescent line of cold iron is one of the most 
powerful probes of cold gas in the circumnuclear envi- 
ronment of an X-ray luminous AGN. The Suzaku spec- 
trum of NGC 4258 shows the most robust evidence to 
date for this iron line in this source (Fig. [2]). Adding 
a Gaussian line to the Canonical Spectral Model leads 
to a very significant improvement in the goodness-of- 
fit (Ax^ = —40 for three additional model parame- 
ters) and gives a line energy, width, flux and equiva- 
lent width of £■ = 6.42 ± O.OSkeV, u < 0.07keV (cor- 
responding to a full width half maximum, FWHM< 
1.1 X lO^kms^i), i^KQ = (6.0li:g) X 10"*^ phcm^^ g-i 
and Wko = 45 dz 17eV, respectively. Assuming Keple- 
rian orbits in an edge-on accretion disk, the limit on the 
FWHM corresponds to r > 3 x 10^ (6 x lO^^pc). 

We note that the equivalent width of the iron line is 
entirely consistent with the limits on reflection reported 
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Fig. 2. — XIS residues from the best-fitting continuum model 
showing the presence of a 6.4 keV fluorescent line of cold iron. 
Data from all XIS are shown (XISO=black, XISl=red, XIS2=green, 
XIS3=blue). 

in Section lOl If, for now, we assume that the iron line is 
produced by isotropic illumination of a planar optically- 
thick structure, we can use the relations reported in 
Matt, Fabian & Reynolds (1997) to infer that the solid 
angle of the reflector as seen by the X-ray source satisfies 
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where 9 is the inclination of the slab, we have assumed 
Anders & Grevesse (1989) abundances. Note that the 
quantity r2/27r can be compared directly with the reflec- 
tion fraction quoted in Section 13.11 A more sophisti- 
cated treatment of the iron line strength expected via 
reflection from the surface of the warped accretion disk 
in NGC 4258 (including all of the relevant geometric ef- 
fects) will be deferred to Section 4.1. 

Once the narrow iron line component has been mod- 
eled, there is no evidence in the XIS spectra for an ad- 
ditional broad/relativistic iron line from the inner accre- 
tion disk. However, due to the high inclination angle 
of the inner disk the limits are not strong. Including a 
Schwarzschild disk-line (modeled using the fully relativis- 
tic code of Brenneman & Reynolds [2006]) with a rest- 
frame energy of -E = 6.4 keV and an emissivity profile 
e ~ between nn = 6rg and Tout = lOOOrg, we derive 
an upper limit on the equivalent width of any broad iron 
line of M4>road < 180 eV assuming an inner disk inclina- 
tion of i = 80°. By contrast, even if the inner disk is 
in an optically-thick fluorescing state and irradiated by 
an isotropic X-ray source, limb-darkening effects would 
likely reduce the equivalent width of the broad iron line 
to below 100 eV (Matt et al. 1997; ROO). Thus, we can- 
not yet rule out the possibility of an optically-thick X-ray 
irradiated inner accretion disk. 

Despite the lower signal-to-noise, the iron line is also 
robustly detected in the new XMM-Newton EPIC data. 
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Fig. 3. — Historical variability of the iron line flux (units 
10~® phcm~^ s~^, equivalent width (units eV), and underlying 6— 
7keV continuum (units 10~^^ erg cm~^ s~^ ). A 10% systematic 
(calibration uncertainty) has been assumed for all continuum flux 
measurements. 

Adding a Gaussian line leads to an improvement in the 
goodness-of-fit of Ax^ = — 26 for three additional model 
parameters, and gives a line energy, width, flux and 
equivalent width of £; = 6.41 ± 0.03 keV, a < 0.07 keV, 
Fkc = (3.3ti;o) X 10-*^phcm-2s-i and Wkc = 53 ± 
19 eV, respectively. Comparing the XMM-Newton and 
Suzaku iron line fits, we see evidence for a decrease in 
the flux of the iron line. Motivated by the desire to use 
line variability to locate the fluorescing matter, this re- 
sult prompts us to conduct a more systematic analysis 
of the long term variability of the iron line. 

3.3. Long term variability of the iron line flux 

Given the weak nature of the iron line in NGC 4258 
(which is itself a comparatively X-ray faint AGN) , there 
are relatively few datasets capable of providing good 
constraints on the line flux. In addition to the deep 
Suzaku and XMM-Newton observations presented here, 
we examine flve additional datasets from the HEASARC 
archives; a deep ASCA observation (15-20 May 1999; 
169 ks of good data), and four shorter XMM-Newton 
observations (8-Dec-2000 [21.6ks]; 6-May-2001 [12.9ks]; 
17-June-2001 [13.6 ks]; 22-May-2001 [16.5 ks]). We note 
that there is also a 15.2 ks XMM-Newton observation 
on 17-Dec-2001 that we choose to ignore due to it be- 
ing severely affected by flares in the instrumental back- 
ground. The XMM-Newton data were processed accord- 
ing to the description given in §2. Processing of the 
ASCA data follows ROO except for the use of the latest 
version of the FTOOLS/HEADAS package (v6.4) and 
calibration files. 

Figure [3] shows the iron line flux, equivalent width, and 
the 6-7 keV continuum flux for all of the data under con- 
sideration. The superior statistics of the new datasets are 
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Fig. 4. — Sum of the 0.3-10 keV light curves from the Front Illu- 
minated XISs (XIS0,2,3). Also shown are the results of a "Bayesian 
block" analysis with p = 3 X 10~^ (orange) and p = 0.085 (red). 



evident. The hypothesis of a constant line flux can be for- 
mally rejected at the 95% level (x^ — 12.1 for 6 degrees 
of freedom), whereas the data are consistent with a con- 
stant equivalent width (a model with Wkcx = 49 eV gives 
= 7.7 for 6 degrees of freedom). The evidence for flux 
variability comes primarily from the new XMM-Newton 
dataset (MJD 54056) which coincides with a historical 
minimum in the 6-7 keV continuum flux. 

Variations of the line flux by almost a factor of two over 
the course of 160 days (i.e., between the new Suzaku and 
XMM-Newton observations) allows us to place interest- 
ing constraints on the location of the fluorescing matter; 
the light crossing time of the full line emitting region 
must be no greater than 160 days (and likely signifl- 
cantly smaller). Thus, we conclude that the line emitting 
material is within a radius of 0.07 pc (4 x lO^r^) from 
the central X-ray source. This strongly suggests that 
there is cold material (producing fluorescent iron emis- 
sion) within the masing portion of the disk (the masers 
appear to have an inner truncation radius of approxi- 
mately 0.13pc). 

3.4. The Suzaku flare 

Significant time variability on a range of characteristic 
timescales is seen during the Suzaku observation (Fig. |4]) . 
A broad lOOksec hump is centered within the (almost) 
200ksec observation window. However, much more rapid 
variability is apparent. Applying a Bayesian Blocks anal- 
ysis^ (based upon earlier work by Scargle 1998), which 
searches for flare-like structure in lightcurves dominated 
by Poisson statistics, the overall lightcurve is found to 
be significantly variable; with robust evidence for 10% 
variations on 16ksec timescales (false alarm probability 
p = 3 X 10~^), and good evidence for the presence of a 
5kscc fiarc {p = 8.5 x 10"2)_ 

^ http:/ /space. mit.edu/CXC/analysis/SITAR/bb_expcriment. html 
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Direct inspection of XIS images confirm that the vari- 
able source is located at the nucleus of NGC 4258 to 
within 20". This confirms and strengthens the case for 
rapid variability seen by BeppoSAX (Fiore et al. 2001), 
and alleviates the background and confusion concerns ex- 
pressed by Fruscione et al. (2005) about the BeppoSAX 
variability. In the absence of relativistic beaming effects, 
light-crossing time arguments allow one to estimate an 
approximate upper limit to the size of the emitting re- 
gions; the large amplitude 50ksec variability should orig- 
inate from a region smaller than about 250rg, whereas 
the corresponding limit for the rapid 5 ksec flare is just 
25rg. If this variability is associated with accretion disk 
processes, more stringent limits can be set by equating 
these timescales with the dynamical time of a Keplerian 
disk = y/r-^/GM. Thus, under this disk-origin sce- 
nario, the 50 ksec variability should originate from within 
r « 40rg whereas the 5 ksec flare should occur in the in- 
nermost disk, r « lOrg. 

We search for any spectral changes during the flare 
by splitting the Suzaku observation into a "low-state" 
and a "high-state" . Only the XIS data have sufficient 
S/N to permit this exercise. We define the threshold 
distinguishing high-state from low state to be a total 
XIS count rate of Icps. Using the Canonical Spectral 
Model, the low-state XISO-3 spectra are fit jointly with 
the high-state XISO-3 spectra. The parameters describ- 
ing the soft thermal plasma emission are fixed between 
the low- and high-state spectral models. Initially, we also 
assume that the photon index of the AGN emission and 
the intrinsic absorption column were also the same for 
the two states, with only the normalizations of the AGN 
power-laws being allowed to fioat between the low- and 
high-state spectral models. The resulting best-fit param- 
eters are identical to that found for the full time-averaged 
spectrum. We then allow the photon index and the ab- 
sorbing column density to be different between the low- 
and high-state. While there was a hint that both the pho- 
ton index and the absorbing column density increased in 
the high-state, the formal improvement in the goodness 
of fit is not statistically significant (Ax^ = —5 for the 
additional of two new model parameters). Thus, direct 
spectral modeling of the low- and high-state spectra does 
not reveal robust evidence of spectral variability. 

An alternative methodology is to construct and then 
model the high— low difference spectrum; this is a for- 
mally correct procedure (and yields meaningful results) 
if the physical difference between the high- and low- 
state is the addition of a new emission component. We 
find that the difference spectrum has the form of a 
pure absorbed power-law; this explicitly demonstrates 
that, as expected, the soft X-ray thermal plasma com- 
ponent has not changed between the high- and low- 
state spectra. Interestingly, both the absorption column 
and photon index characterizing this variable emission 
differ from that found in the time-averaged spectrum, 
iVn = (1.58t[Jj^) X 1023 cm-2 and F = 2.4 ± 0.5 [com- 
pared with A^H = {^■'i.toi) X 10^2 and F = l-l^t^M 
the time-averaged spectrum.] While these results must 
be taken as tentative, it appears that the difference spec- 
tra are revealing spectral changes that are too subtle to 
be revealed by direct modeling of the low- and high-state 
spectra alone. 
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Fig. 5. — Radial distribution of iron line emission in the overall 
best fitting warped disk model of M08 (dashed blue line) and the 
best fitting gravitationally-stable model of M08 (solid red line). 
The quantity shown here /^(r) is the total observed line flux (in 
arbitrary units) from a given radius in the accretion disk. 

4. DISCUSSION 

4.1. The origin of the fluorescent iron line 

Our study has found the first significant evidence of 
iron line fiux variability in this source. Combined with 
the limits on the velocity width of the iron line, the vari- 
ability allows us to constrain the line emitting region to 
the range 3 x lO^rg < r < 4 x lO^'rg (6 x lO'^pc < 
r < 7 X 10~2pc). The most obvious candidate for the 
fluorescing matter on these spatial scales is the accretion 
disk itself. 

Assuming an isotropic X-ray source at the center of the 
accretion disk, the expected iron line from the disk can be 
estimated once we have a model for the 3-dimensional ge- 
ometry of the warped disk relative to the observer. Mar- 
tin (2008; M08) has recently described the maser data for 
NGC 4258 accretion disk using a model of a disk warped 
by the General Relativistic frame-dragging effects of cen- 
tral rotating black hole (Bardeen & Petterson 1975). As 
described in the Appendix, we have calculated the ex- 
pected iron line equivalent width and the radial distribu- 
tion of iron line emission for both the overall best fitting 
warped disk model of M08, and the best fitting warped 
disk model of M08 that is stable to self-gravity. These 
two warped disk models have a total warp angle between 
the inner and outer disk of = 45.8° and 85.9°, respec- 
tively. Assuming solar abundances (Anders & Grevesse 
1989), the = 45.8° model gives an iron line that is 
both too weak (TFkq — 25 eV) and originates from too 
far out in the disk to be compatible with the line vari- 
ability (5 X 10-* - 1 X lO^rg; see Fig. [S]). On the other 
hand, the gravitational stable warped disk of M08 pre- 
dicts an iron line of the correct strength (Wkq = 51 eV), 
originating from a region in the disk that is compatible 
with both the line width and line variability constraints 
((1 - 4) X lO^rg; Fig. [5]). We note that in addition to 
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the high viewing inchnation, the geometry of the warped 
disk imphes that most of the observed iron hne emission 
is driven by highly obhque irradiation of the disk sur- 
face. Thus, it is important to account for the dependence 
of the iron Unc photon production on the angle of inci- 
dence of the irradiating X-ray continuum (see Appendix). 
The high viewing inclination and the oblique irradiation 
means that simple estimates of the line strength based 
on covering fraction arguments are inappropriate. 

Thus, we have shown that one can quantitatively un- 
derstand the properties of the iron line in NGC 4258 if it 
originates from the surface layers of a warped disk irradi- 
ated by a central X-ray source. We do require, however, 
a severe warp inside of the masing region of the disk of 
the kind possessed by the M08 best-fitting gravitation- 
ally stable model. Interestingly, as discussed in M08, 
this model implies a significant mis-alignment between 
the inner accretion disk and the jet. 

The warped disk hypothesis for the origin of the 
iron line in NGC 4258 can be tested by future X-ray 
spectroscopy by searching for the predicted asymmetric 
double-peaked line profile. Fig. [6] shows the predicted 
line profiles from the two warped disk models of M08, 
both for a perfect (5-function emission line at 6.4 keV (left 
panel) and the real Fe Ka doublet with component en- 
ergies of 6.392 keV and 6.405 keV and a 1:2 branching 
ratio (right panel). Future high- resolution spectrographs 
(such as the micro-calorimeters on Astro-H or Interna- 
tional X-ray Observatory) will be able to easily revolve 
the predicted line profiles. 

4.2. The circumnuclear environment of NGC ^258 

As discussed above, most if not all of the iron line emis- 
sion in NGC 4258 can be understood as originating from 
the warped accretion disk. This leaves very little room 
for any iron fluorescence from other (non-disk) cold mat- 
ter in the system such as the putative molecular torus 
of unified Seyfert schemes. Indeed, the circumnuclear 
environment of this black hole does appear to be sig- 
nificantly "cleaner" than the vast majority of Seyfert-2 
galaxies. Dadina (2008) examines a sample of 62 Seyfert 
2 galaxies observed with BeppoSAX and finds a mean re- 
flection fraction oiTZ = 0.87±0.14 (assuming face-on re- 
flection) and a mean iron line equivalent width of Wko = 
693 ± 195 eV. The strong iron line and Compton reflec- 
tion signatures in most Seyfert 2 nuclei are readily under- 
stood as originating from the obscuring molecular torus 
(Krohk, Madau & Zycki 1994). For NGC 4258, the fact 
that our almost edge-on view of the central engine suf- 
fers absorption at "only" the ~ 10^'^ cm~^ level already 
rules out a Compton-thick torus that is aligned with the 
accretion disk. Furthermore, the very weak X-ray repro- 
cessing features {TZ < 0.43 and Wxa = 45 ± 17eV from 
our time-averaged Suzaku data) rules out the presence 
of even a misaligned geometrically-thick Compton-thick 
torus. 

4.3. A comparison with M81* and other AGN 

It is interesting to compare NGC 4258 with the 
LLAGN M81*, since both seem to have very clean nu- 
clear environments. M81 is classifled as a Seyfert 1.8 / 
Low-Ionization Emission Line Region (LINER) galaxy, 
with a central black hole mass of M = 7 x lO^Mo (De- 
vereux et al. 2003) accreting at L ~ lO^^LEdd- M81* 



shows evidence of a one-sided radio jet emanating from 
the nucleus (Bietenholz, Bartel & Rupen 2000). High- 
resolution X-ray spectroscopy of M81* with Ghandra 
(Young et al. 2007) reveals a number of emission lines, in- 
cluding Fe Ka, and velocity broadened Si Ka and other 
thermal fines. The Fe Ka line in M81* has an EW of 
47^24 eV, almost identical to that of NGC 4258, and 
there is no evidence of a broad iron line. Furthermore, 
the broadened Si Ka fluorescence line in M81* is consis- 
tent with originating in a disk at r ~ lO'^r^, assuming the 
inclination angle is the same as the disk observed with 
HST, i = 14° (Devereux et al. 2003). The broadened 
thermal lines in M81* are consistent with originating at 
r < 10''~^rg, suggesting that there is hot thermal gas at 
small radii, possibly in the form of a radiatively ineffi- 
cient accretion flow or the base of a jet (Markoff et al. 
2008). 

The inclination angles of the outer accretion disk in 
NGC 4258 (- 80°), and in M81, (14°) are significantly 
different, and this may account for the fact that NGC 
4258 has a much larger column density (~ 10^^ cm~^) 
than M81 (~ 10^^ cm~^), that maser emission is not 
seen in M81, and the difference in the obscuration of the 
broad line region (i.e., the Seyfert types). In M81*, we 
do not know the geometry of the disk well enough to 
calculate the strength of the Fe Ka line that it would 
produce. Apart from their different inclination angles, 
the accretion flows in NGC 4258 and M81* seem to be 
remarkably similar. 

Although we have noted the similarities between the 
circumnuclear environments of NGC4258 and M81*, it 
is interesting that these results on NGC4258 and M81* 
appear to run counter to some general behaviour of the 
AGN population. In particular, it has been noted that 
the strength of the iron emission line across the AGN 
population is anticorrelated with both the X-ray lumi- 
nosity (Iwasawa & Taniguchi 1993) and the Edding- 
ton ratio (Winter et al. 2008). Using the Swift-BAT 
AGN survey. Winter et al. (2008) show that LLAGN 
with X-ray Eddington ratios comparable to NGC 4258 
{Lx/L-Edd ^ 10^^) typically possess iron line equivalent 
widths in excess of 500 eV. Thus, on the basis of the iron 
line equivalent width, the circumnuclear environments of 
NGC 4258 and M81* appear to have significantly less 
cold gas than the average LLAGN. 

5. CONCLUSIONS 

Using Suzaku, XMM-Newton and Swift, we have ob- 
tained an unprecedented view of the active nucleus in 
NGC 4258. Our principal results are: 

1. Comparing the Suzaku data with XMM-Newton 
data taken 160 days later, we detect robust flux 
variability of the 6.4 keV iron line for the first time. 
Together with constraints on the velocity width of 
the line, and assuming Keplerian motion about the 
central black hole, we can place the iron line emit- 
ting region to be between 3 x lO'^r^ and 4 x lO^r^. 

2. We show that the strength, velocity width and time 
variability of the iron line can be explained by a 
model in which the line originates from the sur- 
face of a warped accretion disk. In particular, we 
present explicit calculations of the expected iron 
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Fig. 6. — Predicted iron line profiles for the best fitting warped disk model of M08 (dashed blue line) as well as the best fitting 
gravitationally stable warped disk model (solid red line). The left panel shows the line profile for a single (5-function emission line at 
(S.4keV; the symmetry of the redshifts and blueshifts is apparent. The right panel shows the expected profile of the Fe-Ka doublet (with 
components at 6.392 keV and 6.405 keV and a 1:2 branching ratio.) 



line from a disk warped by Lens-Thirring preces- 
sion from a severely misaligned central black hole. 

3. During our Suzaku observation, we detect high am- 
plitude intraday variability, with fluctuations on 
timescales as short as 5ks. Corresponding light 
travel time arguments suggest that the emission re- 
gion is smaller than 25rg. If we make the stronger 
assertion that this timescale be longer than the dy- 
namical timescale of the accretion disk at the loca- 
tion it is produced, the upper limit on the radius 
of the emission is lOrg. 

4. In stark contrast with the vast majority of other 
Seyfert 2 galaxies, there are no indications of a 
Compton-thick obscuring torus; the weak iron line 
and the lack of reflection all point to a circumnu- 
clear environment that is remarkable clean of cold 
gas. As pointed out by Herrnstein et al. (2005), 
the intrinsic absorption that we do see in the X- 
ray spectrum may well arise in the outer layers of 
the warped geometrically-thin accretion disk, fur- 
ther reducing the need for any cold structure other 
than the accretion disk itself. 

5. We highlight the similarities in the circumnuclear 
environments of NGC 4258 and another LLAGN, 
M81*. However, we also note that the remarkably 
clean circumnuclear environment found in these 
two LLAGN stand in contrast to the vast majority 
of LLAGN. 

We thank Richard Mushotzky for stimulating conver- 
sations throughout this work. GSR thanks the NASA 
Suzaku and XMM-Newton Guest Observer Programs for 
support under grants NNX06A135G and NNX07AE97G. 
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APPENDIX 

IRON LINE EMISSION FROM A CENTRALLY ILLUMINATED WARPED DISK 

In this Appendix, we describe our calculation of the fluorescent iron line profiles from a centrally illuminated warped 
accretion disk. We base our considerations on the warped disk model for NGC 4258 of Martin (2008; also see Martin, 
Pringle & Tout 2007). Martin et al. (2008; hereafter M08) use the formulism of Scheuer & Feiler (1996) to compute the 
geometry of a steady-state, geometically-thin viscous disk that is warped due to a misalignment between the angular 
momentum of the large scale accretion flow and the spin axis of the central black hole. The interesting parameters of 
the M08 disk model are the misalignment angle between the angular momentum of the large scale accretion flow and the 
spin axis of the central black hole, 77, the warp radius i?warp, and the power-law index describing the radial dependence 
of the surface density b (with surface density following E cx R^^. Setting 6=2, M08 show that such a model describes 
the geometry of the high-velocity masers in NGC 4258. Although some degeneracy exists, M08 described two "best 
fitting" models; the overall best-fit has 77 = 45.8° and i?warp — 1-6 x lO^r^ (corresponding to 5.8 mas), whereas the best- 
fitting model given the constrain of a gravitationally-stable disk is 77 = 85.9° and -Rwarp — 2.47 x lO^rg (corresponding 
to 1.3 mas). 

A basic premise of this approach is that the disk at a given radius is a flat ring centered on the black hole. The 
different annuli exchange angular momentum via a process modeled as shear viscosity. Let l(-R) be a unit vector 
pointing in the direction of the angular momentum of the disk flow at radius R. In a Cartesian coordinate system in 
which the black hole spin is aligned with the z-axis, the fact that the innermost portions of the disk are aligned with the 
black hole spin (Bardeen & Petterson 1975) imply that l(i?) (0, 0, 1) as r ^ 0. We also have l(i?) (sin 77, 0, cost;) 
as r ^ 00, where we have aligned our coordinate system such that the angular momentum of the large scale flow lies 
in the x — z plane. As discussed in Martin, Pringle & Tout (2007) and M08, the disk warp is then described by 



^(^)-r(i) 3V6 [-^) 



(l-i) 



3 \ Rwarp 



(Al) 



where W = + Hy is a complex representation of the unit vector 1, Woo — sin 77, Ki/q is the modified Bessel function 
of order 1/6, and F is the gamma- function. 

Within the context of their models, M08 find good fits to the maser data assuming that both the black hole spin 
axis and the angular momentum of the large scale accretion flow lie in the plane of the sky. Hence, in the Cartesian 
coordinate system defined above, our line of sight to the disk is directed along the y-axis, i.e., in the direction of the 
unit vector e^. 

Given the disk geometry relative to the observer, we now proceed to compute the X-ray irradiation of the disk 
surface and the resulting iron line emission. Suppose that the primary X-ray source is located at the center of the 
accretion disk, and emits continuum photons isotropically at a rate Nc{E)dE — Nq{E / Eq)~^ dE where, from now on, 
F shall refer to the photon index of the power-law source. Following Pringle (1996), we define the angles j{R) and 
/3{R) such that 

1 = (cos 7 sin sin 7 sin /3, cos (3) , (A2) 

i.e., f3 is the local angle of tilt of the disk relative to the z-axis and 7 marks the position of the decending node of the 
disk. We shall also establish a polar coordinate system {R, cj)) on the surface of the warped disk. We shall denote as x 
the position vector of the point on the disk with coordinates (i?, </>); the corresponding unit vector shall be denoted x. 
Pringle (1996) gives 

x = i?(cos 0sin7 -I- sin cos 7 cos /3, sin (/) sin 7 cos (3 — cos cos 7, — sin sin (3) (A-3) 

From eqn. (2.13) of Pringle (1996), the rate at which a patch dS of the disk surface intercepts photons from the 
central primary source (neglecting possible shadowing) is 

'^^^ S^'*'^^' " S^'^'^ ' COS./) sin /3-i?/3' sin 01 (A4) 
where 7 ' = dj/dR and f3' = dj3/dR. The angle between the irradiating flux and local disk normal is given by 

X • dS i?(cos0sin/37 ' — sin(/)/3') 

^""^ " " IdSf " [1 + i?2 (cos sin ^7 ' - sin (j)[3'Y] 1/2 ' ^ ' 
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Given the irradiating continuum photon flux and the incident angle, we use eqns. 4-7 of George & Fabian (1991) 
to compute the number of fluorescent iron line photons generated by the patch, NiRdRd(j). We then apply a scaling 
factor of 1.3 (Matt et al. 1997) to correct the results of George & Fabian (1991) to those expected for Anders & 
Grevesse (1989) solar abundances. Assuming that these photons are emitting isotropically from the optically-thick 
surface of the disk, the iron line photon emission rate from the patch per unit solid angle towards the observer is 
cos9 NiRdRdcp/n, where 9 is the angle between the local disk normal and the observer's line of sight given by. 

• dS — i?(sin(^sin7COs/3 — cos(/iicos7)(7 'cos^sin/3 — /3'sin(/ii) 
IdST" [1 + i?2(cos<^sin/?7 ' - sin(A/3')2]V2 " 

The contribution of a given radial range r — > r + dr to the observed line emission is given by 

MiR) = ^^^Rd^, (A7) 

and the total number of observed iron line photons is Af tot = /o°° -^(^) dR. The observed equivalent width is derived 
by ratioing this against the number of continuum photons emitted per unit solid angle at the line energy, 

Finally, to compute a line profile, we need to know the line-of-sight velocity of a given patch of the disk. Assuming 
Keplerian flow, the velocity of the disk is given by 



v = lxxV-^, (A9) 



and the observed line of sight velocity is 



*^ios = • V = y ^ [sin 4> cos 7 sin^ /? + cos /3(cos0sin7 + sin cos 7 cos /3)]. (AlO) 



